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ABSTRACT  i 

Experiments  were  carried  out  which  subjected  tubes  of  ductile  metals  to  histories  of  straining  in  combined 
tension/compression/torsion  space.  The  computer  controlled  servo-hydraulic  testing  machine  recorded  the 
resulting  stress  histories  together  with  the  imposed  strain  histories.  Intermiuently  along  the  path,  the  straining 
was  halted  to  probe  for  the  current  yield  locus.  The  nose  of  the  yield  locus  in  the  direction  of  straining 
exhibited  quite  accurately  a  circular  geometry,  which  implied  isotropic-kinematic  hardening. 

The  centers  of  the  fitted  yield  circles  in  (a,  'l3  x)  generalized  tension-shear  space,  define  the  back  stress  a 
which  expresses  the  anisotropic  component  of  the  hardening,  a  expresses  the  strength  and  orientation  of  the 
impediment  to  dislocation  flow  through  the  crystal,  i.e.  the  yield  stress  required  to  maintain  plastic  flow. 
Appli^tion  of  a  polygonal  closed  straining  path  caused  the  back  stress  tensor  to  rotate  about  330°  while  the 
material  ro^tion  due  to  shear  was  only  about  1°.  Thus  plastic  spin,  which  in  the  literature  has  been  considered 
to  cause  spin  of  the  plastic  anisotropy,  turns  out  to  be  very  much  smaller. 

Measurements  of  stress  response  to  linear  and  nonlinear  strain  paths  have  been  used  to  generate  the 
evolution  law  for  the  back  stress  and  hence  the  constitutive  relation. 
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A  COMBINED  EXPERIMENTAL  AND  ANALYTICAL  INVESTIGATION  OF  THE 
THEORY  OF  PLASTICITY  INVOLVING  STRAIN  INDUCED  ANISOTROPY 

1.  Introduction 

When  plastic  flow  occurs  in  the  manufacture  of  machine  components  or  structures,  or  is 
anticipated  to  occur  in  the  use  of  them  or  in  associated  strength  and  reliability  design 
calculations,  isotropic  plasticity  theory  has  commonly  been  adopted  because  reliable 
anisotropic  theory  and  the  needed  material  characteristics  are  not  usually  available.  However, 
the  strain-hardening  generated  by  plastic  flow  commonly  produces  anisotropic  yield  surfaces 
and  hence  anisotropic  flow  laws  governing  subsequent  plastic  flow.  The  restriction  to 
isotropic  theory  clearly  leads  to  errors  in  the  prediction  of  the  stresses  generated  by  the  plastic 
straining  imposed  and  hence  prevents  reliable  design  of  the  manufacturing  process  and 
prediction  of  characteristics  of  the  product,  such  as  high  localized  strain  concentrations. 
Experiments  were  therefore  carried  out  on  ductile  metals  at  ambient  temperatures  to 
investigate  the  stresses  and  yield  surfaces  generated  following  imposed  histories  of  strain  in 
tension/compression-torsion  space.  These  are  needed  to  develop  appropriate  constitutive 
relations.  In  order  to  generalize  the  results  to  three-dimensional  theory,  an  additional 
measurement  of  the  varying  diameter  of  the  tube  to  give  the  history  of  circumferential  strain 
was  also  necessary. 

The  experiments  were  carried  out  by  imposing  prescribed  linear  and  nonlinear 
tension/compression-torsion  strain  paths  on  tubular  specimens  using  a  servo-controlled  MTS 
testing  machine.  The  straining  was  interrupted  intermittently  to  probe  for  the  yield  loci  in 
tension/compression-shear  stress-space.  The  whole  process  was  controlled  by  a  computer 
without  the  need  for  human  intervention.  Similarly,  the  vast  amount  of  recorded  data,  such  as 
the  points  delineating  the  yield  loci,  were  plotted  automatically  from  the  computer  output  and 
interpreted  in  terms  of  constitutive  theory  by  the  computer  without  the  intervention  of  human 
judgement. 
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2.  The  Testing  Program  for  Material  Characteristics 

Pure  aluminum,  aluminum-magnesium  and  aluminum-manganese  alloys  were  adopted  as 
examples  of  ductile  metals  which  exhibit  elastic-plastic,  essentially  rate-independent,  material 
behavior,  conforming  with  classical,  strain  hardening,  elastic-plastic  theory.  The  response  to  a 
strain  controlled  tension  test  with  intermittent  unloading  is  shown  in  Figure  1.  It  is  seen  to 
express  elastic-plastic  behavior  with  an  extensive  range  of  strain  hardening.  It  exhibits  an 
intergranular  instability  process  known  as  the  Portevin-le  Chattalier  Effect,  but  with  strain 
control  this  does  not  contribute  a  sufficient  distuihance  to  prevent  application  of  classical 
elastic-plastic  theory.  The  unloading  and  re-loading  exhibit  reversible  elastic  behavior. 
Figure  2  illustrates  tension  test  results  under  stress  control  in  which  the  instability  is  much 
more  marked,  because,  when  plastic  flow  is  intermittently  re-initiated,  the  strain  rate  must 
jump  to  very  large  values  instantaneously  because  the  prescribed  stress  rate  remains  constant, 
and  as  the  tangent  modulus  falls  to  the  small  plastic  modulus  the  plastic  strain  rate  must 
increase  suddenly  which  cannot  be  precisely  accommodated  by  the  servo-drive.  In  contrast, 
under  strain  control  the  very  rapid  reduction  to  the  plastic  modulus  demands  only  that  the 
stress  remains  nearly  constant  and  so  the  stress  rate  drops  suddenly  which  poses  no  problem 
for  the  servo-drive  and  the  material  response.  Strain  control  was  used  throughout  the  testing 
program  in  order  to  avoid  these  instabilities,  and  this  is  in  conformity  with  the  development  of 
elastic-plastic  constitutive  relations  in  which  the  stress  history  is  expressed  in  terms  of  the 
history  of  deformation. 

In  order  to  interpret  tension-torsion  testing,  it  is  convenient  to  utilize  the  generalized 
shear  stress  V3  x,  where  x  is  the  shear  stress  and  correspondingly  the  generalized  plastic  shear 
strain  yP  /  V3.  The  Mises  J2  yield  condition  then  becomes 

CT^+ 3x^=  y^  =  o^  +  (  V3  t  (2.1) 

and  the  corresponding  plastic  work  rate 


G  sF  +  Ty^  =  0^  -»-V3T(y^  /  V3) . 


(2.2) 


The  Mises  yield  condition  is  thus  a  circle  in  (  a,  V3  t )  space  and  the  corresponding 
maximum  plastic  work  concept  demands  normality  of  the  strain  rate  vector  /  V3)to 

the  circular  yield  locus  in  (  a  ,  VSx  )  space. 

The  first  tension/compression-torsion  straining  test  carried  out  followed  a  polygonal  path 
in  ( e  ,  Y  /  V3  )  space.  Probes  for  yield  loci  were  carried  out  at  each  vertex.  Figure  3 
illustrates  the  procedure  for  probing  for  the  yield  locus  of  the  material  at  L,-,  the  vertex  or 
loading  point  station.  The  straining  path  has  been  linear  from  the  station  L,-  _  j.  On  reaching 
Li,  the  straining  is  reversed  and  elastic  recovery  takes  place  to  the  base  point  S,-  inside  the 
current  yield  locus.  A  sequence  of  sixteen  excursions  from  5,-  are  programmed  in  the 
numerical  order  shown  in  Figure  3  to  determine  the  stresses  and  strains  at  the  onset  of 
renewed  plasticity,  that  is  at  points  on  the  new  yield  locus  following  deformation  along  L,-  _  j 
—  Along  each  excursion  the  controlling  computer  records  the  strains  and  stresses  at 
fourteen  points  planned  to  remain  within  the  current  elastic  region  and  carries  out  a  least- 
squares  fit  for  straight  lines  through  both  axial  stress-axial  strain  and  shear  stress-shear  strain 
values  recorded  in  the  computer.  The  plot  for  the  axial  components  of  the  eleventh  excursion 
are  shown  in  Figure  4.  The  linear  excursion  continues  in  strain  space  and  the  computer 
compares  the  stress-strain  evaluation  for  both  axial  and  shear  components  with  the  linear 
extension  of  the  fitted  elastic  lines.  The  strain  deviations  thus  revealed,  Ae  and  Ay  /  V3, 
express  the  renewed  plastic  flow  being  initiated.  An  offset  of  lO"^  was  selected  to  define  that 
a  point  on  the  new  yield  locus  had  been  reached  according  to 

V  (Ae)2  (Ay  /  =  IQ-^  .  (2.3) 

This  offset  magnitude  was  chosen  to  avoid  adding  appreciable  strain  hardening  from  the 
probing  procedure  while  maintaining  a  sufficiently  definite  signal  of  renewed  plastic  flow  to 
define  a  reliable  yield  locus.  For  material  in  the  state  associated  with  loading  to  the  strain 
point  Li,  the  probing  was  repeated  for  each  of  the  sixteen  excursions  indicated  in  Figure  3 
before  continuing  the  program  of  plastic  straining  to  the  next  loading  station.  Figure  5  shows 
the  total  strain  history  to  which  this  specimen  was  subjected.  Some  station  numbers  are 
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indicated.  The  initial  deformation  from  station  1  to  station  2  was  simple  tension  followed  by 
a  sixteen  sided  regular  polygonal  path  with  yield  locus  probing  at  each  vertex.  Finally  a 
square  path  was  followed.  In  some  cases  some  station  strain  locations  were  repeated,  for 
example  2  and  18.  The  end  points  of  the  probing  excursions  delineate  the  yield  loci  in  strain 
space  when  only  one  yield  locus  probing  occulted  at  that  strain,  which  illustrates  the 
extensive  strain  hardening  which  was  generated,  and  this  applies  also  in  stress  space  because 
the  elastic  constants  are  not  changed  by  plastic  flow. 

Figures  6a  and  6b  give  details  of  the  yield  surfaces  measurements  for  stations  6  and  14, 
respectively.  The  broken  lines  indicate  the  direction  of  the  strain-rate  vector  (  e  ,  y  /  V3  ) 
along  the  side  of  the  polygon  traversed  just  before  the  probing.  The  stress  points  shown  are 
automatically  plotted  from  the  computer  output.  From  the  probing  measurements  it  is 
generally  observed  that  the  yield  locus  exhibits  a  protruding  rounded  nose  in  the  direction  of 
the  strain-rate  vector  with  a  flattened  rear  which  involves  appreciable  scatter.  Why  the  latter 
occurs  is  clarified  in  Figure  7.  When  the  probing  is  in  the  general  direction  of  the 
immediately  preceding  plastic  strain-rate  vector,  the  onset  of  plasticity  is  closer  to  that  for 
intermittent  loading  in  a  tension  test,  whereas  for  the  flattened  rear  of  the  yield  surface  the 
situation  is  closer  to  tension  followed  by  compression.  For  a  given  strain  offset  the  Aa  is 
larger  and  more  variable  in  the  case  of  reversed  straining  because  of  the  larger  and  more 
rapidly  changing  tangent  modulus. 

Examination  of  the  yield  loci  established  that  the  measured  points  on  the  rounded  noses 
closely  approximate  a  circular  geometry.  To  check  this  analytically,  least-squares  fits  of 
circles  through  6,  8  and  12  measured  points  on  the  nose  were  evaluated.  Examples  of  such 
circles  arc  shown  in  Figiues  6a  and  6b  and  indicate  a  good  representation  of  the  measured 
points.  Such  circular  yield  loci  constitute  combined  isotropic-kinematic  hardening,  the  radius 
of  the  circle  expressing  the  isotropic  component  and  the  coordinates  of  the  center, 
( o„  ,  V3  ),  the  back-stress  vector  in  tension-torsion  space,  expressing  the  anisotropic 
component  of  the  hardening.  In  general  straining  the  yield  locus  is  a  circle  in  stress  deviator 
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space  with  center  the  back  stress  tensor  a. 

Many  technological  metal  forming  problems  are  governed  by  plastic  deformation 
occurring  on  the  nose  of  the  associated  yield  surface  because  the  objective  of  forming  is  to 
change  the  shape  of  a  workpiece  to  produce  a  prescribed  product  so  that  the  process  is 
continuously  generating  similar  changes  of  geometry.  Thus  the  strain-rate  tensor  is  oriented 
in  a  restricted  sector  of  strain-rate  space  which  forms  the  rounded  nose  of  the  yield  locus.  In 
this  circumstance,  analytical  investigation  of  the  development  of  metal-forming  defects  which 
may  be  generated  during  forming,  can  be  carried  out  on  the  basis  of  combined  isotropic- 
kinematic  hardening.  This  is  convenient  both  from  the  standpoint  of  measuring  the  needed 
physical  characteristics  and  also  of  applying  the  associated  theory. 

Figure  8  shows  the  evolution  of  the  back-stress  vector,  (  Oa  .  )»  centers  of 

least-squares  fitted  yield  loci,  based  respectively  on  the  8  and  12  probed  yield  points  on  the 
rounded  noses  of  the  yield  loci.  The  plotted  points  correspond  to  the  stations  shown  in  Figure 
5,  so  that  the  initial  tensile  straining  and  traverse  of  the  regular-polygonal  path  are  indicated. 
The  initial  yield  locus  at  station  1  shows  that  the  material  supplied  already  exhibited  a  slightly 
anisotropic  yielding  behavior.  The  shear  strain  magnitudes  were  of  the  order  2% 
corresponding  to  a  maximum  material  line  rotation  of  about  1°,  whereas  the  back  stress  a 
rotated  through  about  330°.  Thus  the  rotation  of  the  anisotropic  yield  loci  characteristics 
completely  overwhelms  the  material  rotation. 

3.  The  Structure  of  the  Constitutive  Relation 

The  general  theory  of  elastic-plastic  deformation  at  finite  strain,  with  the  strain¬ 
hardening-induced  anisotropy  modeled  as  combined  isotropic-kinematic  hardening,  was 
presented  in  Agah-Tehrani  et  al.  (J.  Mech.  Phys.  of  Solids,  Vol.35,  pp.5 19-539,  1987).  As 
explained  there,  for  general  plasticity  analysis  it  is  appropriate  to  adopt  the  Kirchhoff  stress  T 
(  detV*  times  the  Cauchy  or  true  stress)  where  V®  is  the  symmetric  elastic  deformation 
gradient.  Since  in  almost  all  problems,  other  than  high-pressure  shock-wave  analysis,  detV^, 
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the  volume  ratio  of  the  elastically  deformed  material  to  the  uns  jessed  material  is  almost 
unity,  T  is  effectively  equal  to  the  tme-stiess  tensor  5. 

The  change  in  the  yield  surface  due  directly  to  plastic  flow  is  termed  strain  hardening. 
It  is  part  of  a  flow  or  incremental  type  process  in  that  the  rate  of  change  of  the  yield  surface 
relative  to  axes  fixed  in  space  is  determined  by  the  rate  of  plastic  straining  plus  a  contribution 
arising  firom  the  spin  of  the  yield  surface  because  it  is  embedded  in  the  deforming  material. 
For  isotropic -kinematic  hardening,  the  back  stress  a  embodies  the  anisotropic  part  of  the 
changing  yield  condition  which  shifts  the  center  of  the  yield  surface  away  from  the  origin  in 
stress  space.  The  size  of  the  yield  surface  expresses  the  isotropic  component  of  the 
hardening.  Thus  the  evolution  equation  governing  the  anisotropic  component  of  the  hardening 
due  to  plastic  flow  must  involve  the  variable  a,  the  total  rate  change  of  a,  less  the 
contribution  of  the  material  spin  W ,  giving  the  appropriate  yield-surface  variable: 

a  -  (  Wa  -  aW^  )  =  d  -  VTa  +  aW' =  a  (3.1) 

which  is  the  Jaumann  derivative.  The  term  (  W'a  -  aW  )  expresses  the  rate  of  change  of  a 
spinning  constant  tensor  (Lee  et  al.,  J.  AppL  Mech.,  Vol.50,  pp.554-560,  1983). 

Plasticity  comprises  the  stress-motivated  flow  of  dislocations  and  other  mechanisms 
through  a  crystal  lattice  which  produces  permanent  deformation.  Strain  hardening  expresses 
the  impediment  to  such  migration  built  up  by  previous  plastic  flow  which  disturbs  the  initial 
level  of  uniformity  of  the  crystal  lattice.  The  anisotropy  of  the  hardening  arises  because  pile- 
ups  of  dislocations  against  barriers  in  the  crystal  structure  impede  flow  more  strongly  in  the 
direction  of  the  flow  which  produced  them.  Thus  the  evolution  equation  for  the  back  stress  a 
expresses  its  rate  of  growth  in  terms  of  the  plastic  strain  rate  and  the  current  value  of  a 
itself.  DP  is  the  macroscopic  quantity  which  incorporates  the  amplitude  of  the  flow  and  its 
direction  while  a  expresses  the  strength  and  orientation  of  the  resistance  to  continued  plastic 
flow.  Since  the  Jaumann  rate  of  a,  the  appropriate  variable  for  the  rate  of  change  of  a,  is  a 
symmetric  second  order  tensor,  it  must  be  representable  by  an  isotropic,  symmetric  tensor 
function  of  a  and  DP .  Thus  the  general  form  which  the  evolution  equation  for  a  must  take 


-  7  - 


is: 

a=h(a  ,DP  )  =  r[^DP  +  a  +  1130^ (a^  -  (tra^)/ /3) 

+  r[4jiaDP  +  D^a  -  2tr(ocD^)//3) 

+  DP  +  DP  a^-  2a(a}DP/3)I )  (3.2) 

where  DP  =  (2tr(D^D^)/3)‘'^  is  the  generalized  scalar  plastic  strain  rate  and  the  Tj’s  are 
functions  of  the  invariants  of  a.  It  is  linear  in  to  preserve  the  rate  independent  plastic 
response  under  consideration.  Special  cases  are  commonly  utilized  in  which  some  of  the  t|’s 
are  taken  to  be  zero.  If  ilj  is  the  only  nonzero  coefficient,  (3.2)  reduces  to  the  rule  first 
suggested  by  Prager  when  he  introduced  the  concept  of  kinematic  hardening. 

4.  Determination  of  the  Constitutive  Relation 

The  plan  of  our  research  program  is  to  carry  out  experiments  in  which  tubes  are 
deformed  in  tension-torsion  space,  involving  a  series  of  prescribed  straining  paths,  and  to 
record  the  resulting  stress  and  yield  loci  histories  in  order  to  determine  the  structure  of  the 
evolution  equation  for  the  back  stress.  Usually  only  one  of  the  terms  have  been  included  in 
applications  of  (3.2).  In  view  of  the  detailed  information  concerning  the  growth  of  the  back 
stress  achieved  on  the  current  project,  the  inclusion  of  three  terms  in  the  evolution  equation  is 
investigated.  For  the  tension-torsion  of  a  thin  walled  tube,  the  evolution  equations  for  the 
back  stresses  (3.2)  can  then  be  reduced  to 

dj  =r\iD^  +  tti  -I-  (ja?  -  jo^  -  12  (41) 

02  =  TiiD?  112^^  02  +  TI3D'’  (- ja?  +  jal  -  12  (4.2) 

“12  =  'ni^?2  +  112^''’ 0112  +  Ti3D'^(ai  -I-  02)012  +  (02  -  Oi)D  12  (4.3) 

when  114  TI5  are  set  to  zero.  Here  subscripts  1,  2,  3  and  12  represent  the  axial, 
circumferential,  radial  and  shear  components  of  o,  d,  DP  and  D  and  D  is  the  rate  of 
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deformation  tensor.  The  experimental  record  consists  of  information  on  the  true  stresses  5^, 
5i2  and  the  logarithmic  strains  e^,  £2  and  £12-  From  these  data,  the  accumulated  plastic  strains 
and  the  accumulated  generalized  plastic  strain  histories  can  be  evaluated.  In  the  (JTj  V3?Ti2) 
stress  space  the  center  of  a  circular  yield  surface  corresponds  to  the  back  stress  components 
(3/2)ai  and  ■^ot22-  Thus,  we  may  curve-fit  a  circular  geometry  through  the  rounded-nose 
portion  of  the  experimentally  determined  yield  surface  and  obtain  and  ai2  from  the  center 
of  the  circle.  All  of  these  results  allow  us  to  tabulate  £f,  £f ,  £f2,  £12,  cxj  and  ai2  against 
at  the  yield  surface  probing  stations.  Here  we  have  employed  the  definitions  =  j  , 

e  =  j  D  (t)dt  and  ^  =  j  {t)dt .  These  are  introduced  as  convenient  mathematical  variables 
in  terms  of  which  to  express  this  incremental  or  rate  type  theory  in  which  the  velocity 
gradient,  D,  is  the  strain  rate  type  variable  used.  They  are  not  finite  strain  components 
because  the  material  rotates  and  the  strain  increments  apply  to  varying  material  directions  in 
the  body.  We  next  smooth  these  tabulated  results  by  fitting  polynomials  through  these  data  to 
obtain  continuous  representations  of  these  quantities  as  functions  of^.  To  extract  the  "n’s  we 
may  recast  (4.1)  to  (4.3)  as 


dai 

d^ 

=  ^11 

d£f 

2 

■  +  Tl2ai  +  Ti3(  ja?  - 

1  2 
jal- 

1  2 
JO-I 

+ 

yaf.)  + 

20i2 

t^l2 

dE^ 

(4.4) 

da2 

d£f 

1  2 

2  2 

1  2 

1  2. 

2012 

d£i2 

(4.5) 

dE^ 

=  Tll 

d^ 

+  1)202  +  Tl3(- jOf  + 

-1- 

yOi2)  - 

d^ 

dttio 

def2 

d£i2 

d^ 

-—Til  + 'n2®12  “nsv®!  +  ®2)ctl2  + 

de 

(02  -  Oi) 

dgP 

(4.6) 

Since  a  is  a  deviator,  can  be  expressed  as 


a3  =  -  (  ttj  +  02  ) 


(4.7) 


All  the  quantities  in  these  three  evolution  equations,  except  the  Ti’s,  02  and  (da2)/(d£^),  are 
now  available.  Since  02  cannot  be  determined,  and  hence  (da2)/(dE'’ ),  from  direct 
experimental  measurements,  we  employ  normality  cast  in  the  form: 
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(d£f)/(dgP)  ^  (d£f2)/(dgP)  (d£|)/(dgP) 

(2/3)  Si  -  tti  Si2  -  ai2  -  (1/3)  51-02 

to  express  02  alternatively  as: 

1  ^  2 

or 

02  =  -  i^i  -  —  (  5i2  -  012  )  (4.10) 

3  dEf2 

This  allows  us  to  infer  an  averaged  value  of  O2  as  a  function  of  ^  from  these  two 
expressions.  To  gain  confidence  in  this  approximate  procedure,  we  estimate  the  conformance 
of  the  associated  flow  rule  from  the  first  equality  of  (4.8).  If  the  associated  flow  rule  can  be 
established  to  within  experimental  scatter,  we  perform  polynomial  fitting  on  the  averaged  (X2 
and  determine  its  slope,  (da2)/(d^),  from  such  a  continuous  representation.  With  this  last 
piece  of  information  determined,  we  can  estimate  the  Tl’s  as  functions  of  if  from  (4.4)  to 
(4.6).  Since  we  also  can  obtain  the  invariants  of  a  as  functions  of  ,  we  can  eliminate  f  in 
order  to  extract  the  structiue  of  the  T|’s  in  terms  of  the  invariants  of  a. 

The  polygonal  straining  path  employed  in  the  experiments  described  above  involves 
comers.  As  the  straining  direction  is  changed  from  one  side  to  the  next  as  a  comer  is 
traversed,  the  direction  of  the  stress  rate  vector  lags  behind  that  of  the  strain  rate  vector  due  to 
the  influence  of  elasticity.  Normality  is  established  only  after  the  stress  rate  vector  catches  up 
with  the  strain  rate  vector.  For  the  case  of  the  polygonal  straining  path,  the  strain  path  length 
of  each  straight  segment  is  not  long  enough  for  this  to  occur.  However,  as  explained  above, 
normality  has  to  be  relied  upon  in  order  to  infer  approximately  the  circumferential  back  stress 
02  due  to  the  unavailability  of  appropriate  experimental  measurement.  Thus  the  measurement 
from  the  polygonal  straining  path  is  not  ideal  for  determining  the  stmcture  of  the  ii’s.  As  a 
result  of  this,  we  selected  another  set  of  straining  paths  which  were  designed  to  minimize 


-  10- 


such  comer  effects.  The  four  straining  paths  in  the  axial-generalized-shear  strain  space, 
consisting  of  two  proportional  straining  paths  and  two  slightly  nonproportional  straining  paths, 
are  shown  in  Figure  9a-d.  The  specimens  were  prepared  from  Al-Mn  and  the  specimen 
designation  associated  with  each  straining  path  is  also  indicated.  The  corresponding  stress 
responses  from  these  four  tests  in  the  axial-generalized-shear  stress  space  are  shown  in 
Figures  lOa-d.  The  measurements  from  these  tests  were  reduced  to  a  form  suitable  for  the 
extraction  of  the  "n’s  as  described  above.  Some  typical  polynominal  fits  of  the  data  are  shown 
in  Figures  lla-b.  These  are  for  the  stresses,  and  of  the  proportional  path  34g2.  Figure 
12  shows  the  trajectories  of  the  centers  of  the  fitted  circular  yield  loci  in  the  axial- 
generalized-shear  stress  space. 

For  each  straining  history,  we  can  solve  for  the  values  of  Tji,  T|2,  and  113  from  (4.4)  to 
(4.6)  at  each  level  of  ^ ,  using  the  continuous  representations  of  the  variables  obtained  from 
the  experimental  data.  The  values  of  “n’s  as  a  function  of  the  invariant  H,  defined  as 
a  2  (ttxx^)‘‘^,  for  these  paths  are  inconsistent.  Attenjpts  to  extract  r]j  and  112  from  the 
experimental  data  by  setting  1)3  in  the  evolution  equation  to  zero  also  did  not  yield  consistent 
results.  We  therefore  solved  for  Tii  from  the  experimental  data  by  setting  both  T|2  and  T|3  in 
the  evolution  equation  to  zero.  The  values  of  were  obtained  from  the  reduced  form  of 
(4.4)  to  (4.6)  for  each  level  of  ^  using  the  least-squares  method.  Figure  13  shows  the  plot  of 
Tjj  versus  H  for  the  four  straining  paths.  The  four  curves  do  not  collapse  into  one.  However, 
the  two  proportional  straining  paths,  34g2  and  34g3,  do  show  a  certain  degree  of 
correspondence.  Thus  we  curve-fitted  the  values  of  Tii  from  these  two  sets  of  proportional 
straining  data  to  obtain  the  relation 

'nj(  H  )  =  mg  +  n  m2 

where  mg  =  1104.2MPa,  mi=-11.3  and  m3  =  -  0.15MPa“^  Using  such  a  function 
representation  for  Tij,  we  went  back  to  the  data  from  the  nonproportional  straining  paths  34gl 
and  34f3  and  solved  for  r\2  and  113  from  (4.4)  -  (4.6)  using  the  least-squares  method.  The 
resulting  r\2  and  T|3  obtained  from  the  340-straining-path  data  are 
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Tl2(  H  )  =  no  ®  +  ”2 

'n3(  H  )  =  ^0  +  ^  1  ®  +  ^2 

where  no  =  76.1MPa,  ni=-2A,  ^2  =  "  0-001 9MPa-^  ito  =  7.16MPa,  iti=-0.48  and 
^2  =  0.0092MPa"^  The  data  from  the  other  nonproportional  straining  path,  34gl,  did  not 
yield  a  satisfactory  result.  Inspection  of  the  stress  history  of  straining  path  34gl  in  Figure  10a 
shows  a  sharp  bend  in  the  curve  immediately  following  the  elastic  response.  We  cannot  trace 
the  cause  of  such  an  unexpected  behavior.  Hence  the  data  from  this  nonproportional  straining 
path  might  not  be  reliable. 

There  are  thus  two  sources  of  complexity  involved  in  our  attempt  to  extract  the  structure 
of  the  functions  Tii,  TI2  and  '113.  Firstly,  the  circumferential  back  stress  had  to  be  inferred  from 
normality  rather  than  determined  from  direct  measurements.  Secondly,  inspection  of  the 
recorded  stress  histories  in  Figures  lOa-d  indicates  that  the  response  started  to  become  jerky 
midway  through  each  straining  history.  In  future  work,  it  would  be  advisable  to  seek  the 
possibility  of  obtaining  more  complete  experimental  data  so  that  the  use  of  normality  can  be 
avoided.  It  would  also  be  advantageous  to  reduce  the  noise  level  of  the  stress  response  in  the 
experiments. 

5.  Conclusions 

The  findings  of  this  project  have  established  that,  for  some  ductile  aluminum  alloys,  the 
very  significant  anisotropy  of  the  plastic -flow  characteristics  which  develop  during  plastic 
straining  can  be  quite  accurately  expressed  by  combined-isotropic-kinematic  hardening  theory. 
This  applies  for  continuing  plastic  flow,  such  as  occurs  in  many  metal-forming  processes. 
The  continuum  mechanics  theory  of  combined-isotropic-kinematic  hardening  is  well  developed 
for  stress  and  deformation  analysis,  which  thus  offers  the  opportunity  to  include  the  study  of 
the  influence  of  anisotropic  material  response  in  forming  analyses. 
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The  adoption  of  tension/compression-torsion  testing  measurements  provided  the 
opportunity  to  include  nonproportional  two-dimensional  straining  trajectories  to  provide  input 
for  developing  appropriate  constitutive  relations.  Since  plastic  deformation  is  essentially 
incompressible,  the  material  tensors  which  occur,  such  as  the  back  stress,  are  deviatoric  and 
so  can  be  adequately  investigated  using  two-dimensional  testing.  However,  in  order  to  deduce 
the  appropriate  constitutive  theory  it  is  essential  to  carry  out  such  measurements.  In  the 
engineering  literature  it  has  been  common  to  measure  the  yield  stresses  in  tension  and 
compression,  and  if  they  differ,  simply  adopt  an  isotropic-kinematic  hardening  law  determined 
by  these  two  yield  points.  The  present  study  brings  out  the  importance  of  more  generally 
based  two-dimensional  testing  in  order  to  deduce  reliable  constitutive  relations  for 
applications. 


CONTROLLED  TENSION  TEST 
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PROBING  FROM  BASE  POINT  WITH  OFFSET  V(  Ae (  Ay  /  V3  )2  >  10"  ^ 


Figure  4  An  example  of  the  yield  locus  probing  procedure 
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AXIAL  STRAIN  (%) 


Figure  5  The  prescribed  strain  history  of  the  polygonal  and  square  paths 
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Figure  9a  Prescribed  nonproportional  straining  path  34gl 
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(%)  (e)x^^)s  /  Niv>ij.s  wans 


Figure  9b  Prescribed  proportional  straining  path  34g2 
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Figure  9c  Prescribed  proportional  straining  path  34g3 
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Figure  9d  Prescribed  nonproportional  straining  path  340 
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Figure  10a  Recorded  stress  response  from  nonproportional  straining  path  34gl 
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Recorded  stress  response  from  proportional  straining  path  34g3 
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Recorded  stress  response  from  nonproportional  straining  path  34f3 
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Figure  lib  Shear  stress  versus  generalized  plastic  strain  from  34g2 
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